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We study the spontaneous emission of a partially excited 
Bose-Einstein condensate composed of two- level atoms. The 
formation of polaritons induced by the ground-state part of 
the condensate leads to an avoided crossing in the photon 
spectrum. This avoided crossing acts similarly to a photonic 
band gap and modifies the spontaneous emission rate. 



It is well-known that the radiation properties of atoms 
can dramatically be manipulated by changing the envi- 
roment where the atoms emit photons. For microcavities 
it has been demonstrated and for periodic dielectric 
media predicted [|]|| that the supression of spontaneous 
emission (SE) can be quite high. In the case of a mi- 
crocavity this happens because its geometry reduces the 
radiation-mode density, whereas in a periodic dielectric 
medium SE is supressed due to the formation of photonic 
band gaps (PBG). 

The recent achievement of Bose-Einstein condensation 
in magnetic traps Q has provided a new state of mat- 
ter where all atoms share a single macroscopic quantum 
state. Such a state of matter offers great opportunities to 
explore and test new phenomena related to macroscopic 
quantum coherence. Recently several authors have theo- 
retically studied the optical properties of a Bose-Einstein 
condensate (BEC). For example, in a homogeneous, ex- 
tended BEC polaritons (superpositions of excited atoms 
and photons) are formed ||. In a spatially confined 
BEC, the continuous center-of-mass momentum distri- 
bution leads to an increase of SE In addition, the 
stimulated emission can be increased by the Bose en- 
hancement in a BEC § . In the case of two BECs, inter- 
ference effects can be important . 

In the present paper we focus on the study on how the 
macroscopic quantum coherence of a ground-state ho- 
mogeneous BEC can modify the radiation properties of 
atoms which are partially excited from the BEC by a co- 
herent light field. In difference to a microcavity and PBG 
materials, the manipulation of the radiation properties in 
a homogeneously extended BEC is caused by the forma- 
tion of polaritons due to the existence of ground-state 
macroscopic quantum coherence. 

The model: We consider a BEC composed of two- 
level atoms which is coupled to the electromagnetic field. 
The coupling is described by using minimal coupling 
in rotating-wave approximation under neglection of the 



term quadratic in the electromagnetic field. The interac- 
tion Hamiltonian is then given by 

H int = J d 3 kd 3 k'C a {k)a a (k)^ g (k')^t{k + k') +E.C. , 

(1) 

with Ccr(fc) : = uad ■ e cr (fc)[?i/(2(27r) 3 £oWfe] 1 / 2 for a elec- 
tromagnetic mode with polarization vector e a (k). The 
Heisenberg equations of motion for the photon annihila- 
tion operators a a (k) and the field operators ^f e and $? g 
for excited and ground-state atoms can be derived easily 
and are given by 
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ih^ g (k) = 



h z k 

d s k'J2C(k')aUk')^e(k + k') (3) 



iha a (k) = huj k a a (k) + Uk) / d 3 k'^Jk')^ e (k + k') (4) 



We have neglected the nonlinear interatomic interaction 
terms as they are small compared to the coupling energy 



hy/ujQiSg, where we have defined v g := \d\ 2 p g /(2hso) with 
p g being the density of atoms in the ground-state and d 
denoting the atomic dipole moment. 

Macroscopic coherent solution: First, we consider the 
physics in the coherent regime where the atoms in the in- 
ternal ground state share a single macroscopic wavefunc- 
tion ip g and interact with a classical laser field a£. lass (fc). 
This interaction prepares the BEC into a coherent super- 
position of the internal ground-state ip g and an excited 
state ip e . The equations of motion for these three col- 
lective fields are given by Eqs. (||) to (|J) if the quantum 
field operators are replaced by the corresponding classical 
coherent fields. 

We are interested in finding a particular solution of 
these equations corresponding to a BEC coherently cou- 
pled to a running laser wave. We thus make the ansatz 
ip g (k) = (2ir) 3 / 2 ^~pg~5(k) exp[— ipt] for a homogeneous 



ground-state BEC of density p g , a„ s (k) — exp[— iujT,t] 
S(k-k L ) S ai(TL n L [2(2Tr) 3 he uJk L ] 1/2 /{\d\Lu ) correspond- 
ing to a laser beam with frequency lor,, Rabi frequency 
CIl > 0, polarization ctl, and wavevector fcj, {inside the 
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BEC),and<i/>e(fc) = {2i:f /2 ^p~ e 5{k-k L ) exp[-i(/i+w L )t] 
describing coherent excited atoms of density p e and of 
momentum hk^. Inserting these expressions into the 
Heisenberg equations of motions leads to a set of alge- 
braical conditions which fix the parameters p, fcz,, and 
p e which we assume to be smaller than p g . If we ne- 
glect the kinetic energy the density of excited atoms is 
given by «J~pl = ^/p^VLl/{p + A), where A := ujl - 
is the detuning of the laser beam. The wavevector is 
fixed by to kL = uj l /2 + ^/(u L /2) 2 - u^vjjjx + A). For 
A < (A > 0) the chemical potential is given by 
p = -A/2 ± ^(A/2) 2 + Q, 2 L which implies p + A > 
(p+A < 0), respectively. Note that for A > the expres- 
sion for u>k L implies the additional constraint p+A > Av g . 

Derivation of polariton modes: Having found the 
macroscopic coherent solution we now address the ques- 
tion of how quantum fluctuations will affect it. The cor- 
responding stability analysis will determine how the SE 
process occurs as the system evolves away from the ini- 
tial coherent solution. To perform this analysis we lin- 
earise Eqs. (||) to (|J) around the coherent fields by in- 
serting = exp[— ipt]{%jj g + 8^g\ and 9 e — exp[— i(p + 
u)t]{ip e + S^e} as well as a a = exp[— iwi]{a£. lass + 5a a } 
and retaining only terms linear in and 8a a , which 
describe the quantum fluctuations around the coherent 
solution. 

The resulting linearized equations of motions can also 
be derived from an effective Hamiltonian for the quantum 
fluctuations, 



-fffluct — H po \ + -ffspont 



(5) 
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S acT (k)S^l{k)+R.c. } (6) 



describes the formation of polaritons and conserves the 
total number of particles. In contrast to previous work 
||, it also includes an interaction with the ground-state 
fluctuations S^ g mediated through the laser field Ql- 
But since the interaction energy Ti^J VgUj^/LOk, (which we 
from now on will approximate by Ti,Ji>gUJo) is much larger 
than hflo it will only give small corrections to the polari- 
ton energy. We therefore can neglect it. For the same 
reason we will omit the kinetic energy h 2 k 2 /(2M) of the 
particles as well as the chemical potential p. 

The eigenmodes of -ff po i are then easily found. To 
derive them it is of advantage to introduce the "non- 



coupled" polarization vector £*nc(p) which is perpendic- 
ular to the photon momentum p and the atomic dipole 
moment d. Since the interaction couples the polariza- 
tion to d a mode proportional to this vector describes 
uncoupled photons which do not contribute to SE. A sec- 
ond, "coupled", polarization vector £*c(p) is introduced 
by being perpendicular to p and £nc(p)- It is easy to 
show that the scalar product £c(p) ■ d is then given by 
|d|sinz^, where dp is the angle between p and d. We 
associate with £c(p) the quantum fluctuation operator 
Sa(k) := 8a,j = c{k). 

The relevant eigenmodes of -ff po i then consist of two 
parts, (i) free ground-state atoms ip g ^(p) = 5(p — k) with 

momentum hk and energy E g ~ 0, and (ii) polaritons 
with frequency spectrum 



w fc,± = - A ± w k 



(7) 



with Wj: 



yj (A fc /2) 2 + VgiOQ sin 2 and A k 



: Wfc — U>Q. 



The polariton modes are given by 

= S(p- k) f ^/JTguJ^ sin tip 



wp ± Ap/2 



(8) 



where the first component refers to 8^> e (p) and the sec- 
ond one to Sa(p). The polariton spectrum 0J±,p clearly 
exhibits an avoided crossing around Ap = of width 
y/VgU>o sin dp (see Fig. 1). It also contains a small gap 
whose edge is reached in the limit — > and oo || . 
The second part of the Hamiltonian, 



-^spoilt 



d 3 k 



imti k ^5a\k)S^l(k L - k) +H.c.} 



(9) 



describes the spontaneous decay of excited atoms in the 
coherent field. It depends on the density of coherently 
excited atoms p e through v e := p e \d\ 2 /{2he )- Since it 
has a structure similar to the exponent of a squeezing 
operator we can expect that the coherent field configura- 
tion (corresponding to the vaccum plus an external field) 
is instable. 

Spontaneous emission of polaritons: Having derived 
the polariton modes we are in the position to calculate 
the SE rate. This is done by expanding the field operators 
in terms of the eigenmodes, 

(<5* e (fc), Sa(k)) = fd 3 p^2 <PnAk)K,ff ■ (10) 

■* n=± 

The operator b n p annihilates a polariton in the associ- 
ated mode. In this base -ff po i can be written as 



■> n=± 



(11) 
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The interaction hamiltonian H spont which is responsi- 
ble for SE can be written as 



-ffspont = huJoV^ ( —7= ^ sinl9 

J y/uJ^ ^ 

v n=± 

x |6^^|(fci-fc)+H.c.| 



2w k 



(12) 



To derive the SE rate we have to determine the time 
evolution of polaritons around the coherent field config- 
uration. We assume that initially all atoms and photons 
are described by the macroscopic coherent solution, or in 
other words, the quantized polariton field is initially in 
the vacuum |0). This state then evolves under the ac- 
tion of the polariton Hamiltonian. To describe this state 
\ip{t)) we make the following ansatz, which corresponds 
to the one-photon approximation, 



n=± 



(13) 



The Schrodinger equation ih\ip) — Ha uct \ij}} then can 
be solved by using the Laplace transform a(s) — 
J °° exp[— ts]a(t)dt. In a standard manner j|] we then 
arrive at the equation for the Laplace transform, 



a(s)=a(0)[s + I(s)] 1 
I(s) can be written in the form 



(14) 



sin 2 $g 



/ . 1 • 1 v g uJa sin 2 



(15) 



where V denotes the quantization volume. This integral 
agrees with the one found in absence of a BEC (which 
describes SE in free space) by setting v g = A = 0. We 
denote this free space integral by Iq := I(u g = A = 0). 
Both integrals are linearly divergent and can be treated 
in the way pointed out by Bethe (see, e.g., Ref. [flo||), 
i.e., we renormalize the integrals by subtracting the free- 
electron contribution, 



rR.cn 



= I - 



Vv e U. 



(27T) 
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d 3 k- 
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(16) 



and /^ cn 



I Ren (v g = A = 0). These renormalized 
integrals are only logarithmically divergent. 

At this point it is customary in the calculation of the 
free-space SE to perform the Wigner-Weisskopf approx- 
imation by neglecting the dependence of I^ cn (s) on s. 
In the presence of a band gap this is inappropriate due 
to the strong variation of the mode density around the 
gap Nevertheless, we can perform a generalized 



Wigner-Weisskopf approximation in the following way. 
We expect that the typical timescale on which SE hap- 
pens is much larger than the optical cycle timescale l/ui . 
From the definition of the inverse Laplace transform, 



a(t) 



1 

2vri 



6 + 200 



e ts a{s)ds 



(17) 



it is clear that the variable s plays more or less the role 
of a Fourier-transformed time. We thus expect that only 
values of s much smaller than loq contribute significantly 
to the SE. This implies that we can neglect s whereever 
it appears together with luq or u>l- Thus, we are allowed 
to set s — iujl ~ —i^L in the denominator of I(s) while 
retaining the term depending on s — iA. 

In the case of J Ren this procedure immediately repro- 
duces the Wigner-Weisskopf result J^ cn ps AT e {(7o/2) + 
iA 2 "^ b }, where 70 := (Puj 3 / (3irhs c 3 ) denotes the SE 
rate in free space and N e :— Vp e the number of excited 
atoms in the BEC. To fix A 2 "^ b we follow the theory of 
Bethe (see, e.g., |Q) and introduce a cut-off frequency 
of m e c 2 /h in I^ cu , where m e is the electron's mass. Cal- 
culating the principal value of the integral then leads to 
Z^Lamb s» 270. In contrast to free space the SE rate in a 
BEC depends on A 2 "^ b since such a radiative frequency 
correction shifts the center of the avoided crossing (or of 
a band gap ||). 

It remains to calculate a renormalized expression of 
the integral J Ren in the presence of a BEC. Fortunately, 
this task reduces to integrals proportional to I^ cn and 
a couple of convergent integrals. Within the generalized 
Wigner-Weisskopf approximation we find 



rRen I 1 J — « \ rttcn ■ - ■ c /u i - ■ y 

5{s-iA)J 5th \l5(s-zA) 
8 (s - iA) 
3 



rR.cn 



iV e 7o / 47w s 



IV a 



1 _ 4w s 



(s - iA) 



:(1+ (S ' A) ) 3/2 arcoth(J 

iv g y 



1 + ^- } -)\ (18) 



To perform the inverse Laplace transformation we have 
to know the analytical structure of a(s) of Eq. (|TJ). Be- 
cause of the complicated structure of Eq. (^8[) we even- 
tually have to rely on numerical methods to determine 
a(t), but much can be gained by doing a general analy- 
sis first. In general, a(s) has several poles and a branch 
cut originating from the term including the arccoth in 
7 Rcn . This cut lies between s = iA and s = i(A — v g ). 
The inverse Laplace transform a(t) is then given by the 
residues of a{s) at the poles Sj plus a contribution from 
the branch cut, 



a(t) = a 



E 

poles 



e tSi Res(a(s), s — > Sj) 
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lim 



ds 



[a(s + e) - a(s - e)] , (19) 
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where we have used that the poles seem always to be of 
first order and only the straight parts of the integration 
around the branch cut do contribute. We see that to 
each pole corresponds a fraction of excited atoms whose 
decay rate is given by the real part of the pole Si. The 
branch cut corresponds to a fraction of excited atoms 
which decays in a non-exponential way. 

The magnitude of the SE modification is determined 
by the ratio of v g to \s — zA| which in free space is of 
the order of 7o-/V e . Since for contemporary BECs v g is at 
best in the order of 70 SE is notably changed only if there 
are very few excited atoms N e , which implies |A| ^> Vt^. 
As discussed above, the macroscopic coherent solution 
implies in this case for A > the additional constraint 
A > Av g . In Fig. 2 the real part of the two dominat- 
ing poles s±, S2 is shown as a function of A for the case 
N e = 1 and v g = 2.570 (corresponding to an atom den- 
sity p g of 5 x 10 14 cm -3 ). For A > a third pole appears 
with a very small negative decay rate (< 10~ 3 7o). The 
occurence of negative decay rates is consistent within the 
range of validity of the linearized equations for the quan- 
tum fluctuations and may indicate the formation of an 
atom-photon bound state Obviously the change in 
the SE rate can be quite large for small |A|. According 
to Eq. (|l9|) the fraction of atoms belonging to the poles 
can be easily calculated by determining the residue at 
the poles. It turns out that the poles whose real part 
asymptotically approaches 70 always dominate and that 
the fraction of atoms belonging to other poles is signifi- 
cant only for small |A|. The same holds for the fraction 
corresponding to the branch cut. 

If |A| » 70 holds the dominant pole s\ can be cal- 
culated by perturbation theory. Its real part (the decay 
rate) is then given by 



(20) 



We see that the SE rate is altered by a factor of 1 — 
4f g /(5A). It depends on the sign of A whether SE is 
increased or decreased. 

In conlusion we have shown that for N e = O(l) the 
spontaneous emission is suppressed by the presence of a 
BEC in the internal ground state. The latter induces 
an avoided crossing in the polariton spectrum which acts 
similarly to an effective band gap, thereby suppressing 
spontaneous emission. 
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FIG. 1. A BEC induces an avoided crossing in the polariton 
spectrum. Far away from the avoided crossing the polaritons 
describe excited atoms or photons. Thus, if one focuses on 
the photons, the avoided crossing provides an effective band 
gap. 
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FIG. 2. Spontaneous emission rate of a partially excited 
BEC for detuning A > 4f 9 (solid lines) and A < (dashed 
lines). The quasi-particles break up into different fractions 
with different decay rates. The dominating fraction is the 
one whose decay rate asymptotically approaches 70 . 
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